Vacuum solutions for scalar fields confined in cavities 
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We look at vacuum solutions for fields confined in cavities where the boundary conditions can rule 
out constant field configurations, other than the zero field. If the zero field is unstable, symmetry 
breaking can occur to a field configuration of lower energy which is not constant. The stability of the 
zero field is determined by the size of the length scales which characterize the cavity and parameters 
that enter the scalar field potential. There can be a critical length at which an instability of the zero 
field sets in. In addition to looking at the rectangular and spherical cavity in detail, we describe a 
general method which can be used to find approximate analytical solutions when the length scales 
of the cavity are close to the critical value. 
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I. INTRODUCTION 

The idea of spontaneous symmetry breaking plays a 
central role in our current understanding of the standard 
model. (For an excellent treatment by one of the ma- 
jor contributors see Ref. [Q for example.) A toy model 
potential demonstrating the most basic feature is the 
double- well potential for a re al scalar field considered by 
Goldstone f|. (See Eq. (U) below.) This potential has 
minimum away from tp = at non-zero values of tp = ±o, 
and it is easy to show that these non-zero values represent 
possible stable ground states for the theory. 

For flat Minkowski spacetime, the homogeneity and 
isotropy of the spacetime imply that the ground state 
must be constant. However it is easy to envisage a sit- 
uation in which this is no longer the case, even for ho- 
mogeneous and isotropic spacetimes. A simple example 
was provided by Avis and Isham |1 where the spatial 
section of the spacetime was a torus. In this case, since 
a real scalar field is regarded as a cross-section of the 
real-line bundle with the spacetime as the base space, 
there is the possibility of a twisted field. An equivalent, 
but simpler, observation is that we could choose antiperi- 
odic boundary conditions for the field on the torus just 
as well as periodic boundary conditions. Antiperiodic 
boundary conditions mean that the only constant field 
solution consistent with the boundary conditions is the 
zero field. For the double-well potential, this raises the 
interesting question of what the ground state is if cannot 
correspond to a minimum of the potential. What Avis 
and Isham |3| showed was that if the torus was below 
a critical size set by the parameters of the double-well 
potential, tp = was the stable ground state; however, 
if the torus was larger than the critical size, tp = was 
not stable, and the stable ground state was a spatially 
dependent solution. Thus the boundary conditions can 
prohibit constant, non-zero field solutions, and make the 



determination of the ground state more complicated. 

The main purpose of the present paper is to show that 
the situation for scalar fields confined in cavities is similar 
to that studied by Avis and Isham || for twisted scalar 
fields. In Sec. II we will present a fairly general analysis 
showing how to characterise the critical size of cavities 
corresponding to the stability (or instability) of the zero 
field as the ground state. In situations where ip = is 
unstable, an exact determination of the ground state is 
not possible in general. However, when the length scales 
of the cavity are close to the critical values at which an 
instability of <p = sets in, it is possible to find approx- 



imate analytical solutions; this is described in Sec. Ill 



One case where an exact analytical solution is possible 
occurs for a rectangular cavity where the field vanishes on 
only two oppo site p airs of t he ca vity walls. This case is 
studied in Sec. IV A. In Sec. IV B| we use the approximate 
method of Sec. [II to obtain a solution which is shown to 



agree with a n expans ion of the exact solution we found 
in Sec. IV A . In Sec. IV C we study the case where the 
field is totally confined by the cavity, vanishing on all of 
the box walls. The case of a spherical cavity, with the 
field vanishing on the surface of a spherical shell, is ex- 
amined in Sec. [v]. No exact solution could be found, but 
we obtain a numerical solution, as well as apply the ap- 
proximate method of Sec. HI. The final section contains 



a brief discussion of the results. 



II. GENERAL FRAMEWORK 

We will consider a general region of D-dimensional Eu- 
clidean space S with boundary 9S. (The basic formalism 
we will describe does not require the space to be flat, but 
we will not consider curved space in this paper.) The 
action functional for a real scalar field tp will be chosen 
to be 
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SM = J 



dt I da x \ -d^ipdup -U(<p) 



(2.1) 
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where U (<p) is some potential, and da x denotes the in- 
variant volume element for £ which of course depends 
on the choice of coordinates. We will concentrate on the 
case of the double-well potential 



2\2 



(2.2) 



in this paper, although the analysis is easily extended to 
other potentials in a straightfor ward manner. The field 
equation which follows from Eq. (2.1) is 



a<p + U'(<p) = 



(2.3) 



it can be seen that if any of the eigenvalues A n are nega- 
tive, then E[ip + ip n ] < E[ip] and we would conclude that 
the solution tp is unstable to small perturbations since 
a solution of lower energy exists. If all of the eigenval- 
ues are positive, then tp is stable to small perturbations 
(locally stable). 

Su ppo se that we consider the sol ution p = to 
Eq. (2.5) wit h th e potential Eq. ( |2.2| ). The eigenvalue 



equation Eq. (2.8) becomes 



-V 2 - ^a z ) i •„ = A„ i •„ 



(2.9) 



The aim is to solve this for tp in the region £, subject to 
some boundary conditions imposed on tp on the bound- 
ary, t?£. 

The energy functional (which is just the spatial integral 
of the Hamiltonian density) is 



E[<p] = J do x j^ 2 + i|V^| 2 + U{fp) 



(2.4) 



Because the spacetime is static, if we concentrate on the 
lowest energy solution, or ground state, we would expect 
it to be static as well, resulting in no contribution to 
the kinet ic energy term in Eq. (2.4). For static solutions 
Eq. (2.3) becomes 



V 2 <p + U'(<p) = 



(2.5) 



If we assume that (p has no time dependence, then 
Eq. (O) reads 



E[<p] 



)<p + U{<p) 



(2.6) 



after an integration by parts, provided that the boundary 
conditions are such that </?V ' n tp vanishes on <9£, with V„ 
the outwards directed normal derivative, (ie V n = n ■ V 
with n a unit normal vector on <9£ which is directed out 
of the boundary.) We will assume that this is the case 
here. 

For the potential Eq. (2.2), it is obv ious that ip = 
and tp = ±a are solutions to Eq. (2J3). There are two 



key issues that arise at this stage. The first is whether 
or not the boundary conditions on dT, are satisfied; if 
not, then the solution must be rejected. The second is, 
even if the boundary conditions are obeyed, the solution 
must be stable to pert urb ations. Assuming that we have 
a solution ip to Eq. ( |2.5| ) which satisfies the boundary 
conditions on c*X, we can look at the energy of a per- 
turb ed s olution tp + ip where tp is treated as small. From 
Eq. (§j|) we find 



E[tp + ip]-E{ 



da x ip i 



U"i<p))l> (2.7) 



if we work to lowest order in ip. By considering the eigen- 
value problem 



(-V 2 + U"(tp)) ip n = \ n ip n 



(2.8) 



If we consider the lowest eigenvalue Ao, whose sign deter- 
mines the stability of p = 0, this sign requires knowing 
the smallest eigenvalue of the Laplacian —V 2 . Let 



V 2 Vo = fah 



where £q > and is real. We then have 



Ao = i\ ~ \a 2 
o 



(2.10) 



(2.11) 



If the boundary conditions on p (and hence tp n ) are such 
that £q = 0, then we can conclude that Ao < and so 
ip = is unstable. This is the case where the boundary 
conditions allow constant values of the field, and is the 
case for flat Euclidean space. However, if the boundary 
conditions prohibit constant values of the field, we will 

is 
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have £q > 0, and the stability or instability of p 
determined by the magnitude of l\ 



in relation to -^a' 



Typically £q will depend on the length scales describing 
£ and its boundary, and so we will find critical values 
of these lengths at which tp = becomes unstable. The 
ground state will not be constant in this case. The first 
demonstration of this was given in Ref. Q for the case 
where £ was a circle, or a torus, with the field ip satisfying 
antiperiodic boundary conditions. We will see examples 
of this for fields in cavities later in this pap er. 

The other constant solutio ns to Eq. (2__5) are p — ±a. 
The eigenvalue equation Eq. (2.8) reads 



-V 2 + ^a 2 ) ip n = X n 1p n 



(2.12) 



for these solutions. In place of Eq. (2.11) we have 



An = 



A 2 
—a 

3 



(2.13) 



with £q defined as in Eq. (2.1C). Provided that the bound- 
ary conditions allow constant values of the field, if tp = 
is unstable we see that Ao > implying stability of the so- 
lutions tp — ±a. This is the usual situation in Euclidean 
space Q. However if the boundary conditions rule out 
p = ±a as possible solutions, then the ground state when 
tp = is unstable necessarily involves non-constant fields. 

A simple example where constant field values (other 
than zero) are not allowed occurs for a field confined to a 



3 



cavity E with Dirichlet boundary conditions imposed on 
the boundary dT, of the cavity, [ie (p = on 9E.) In this 
case, if <p = is not stable we are faced with the question 
of what the ground state is. In Sec. we will look at the 
simplest case of a cavity which is rectangular. In Sec. |v| 
we will study what happens for a spherical cavity. The 
case of the rectangular box allows for an exact solution for 
the ground state for some boundary conditions; whereas 
for the spherical cavity we were not able to solve for the 
ground state except numerically. Before proceeding to 
these examples, we will present a method for obtaining 
approximate analytical solutions for the ground state. 



III. APPROXIMATE ANALYTICAL RESULTS 

In Ref. [[l| a systematic method was described for ob- 
taining approximate solutions for the ground state in 
cases where it was not possible to solve Eq. ( |2.5| ) exactly 
(except numerically). In this section we will present a 
variation on this method applicable to general cavities. 
We saw that the stability of <p — relied on the mag- 
nitude of £q which in turn depends on the length scales 
of the cavity. For simplicity assume that there is only 
one length scale L in the problem. Let L c be the critical 
value of this length defined by the condition 



A 2 

—a 

6 



(3.1) 



On dimensional grounds, we must have i\ cx L~ 2 , so for 
L < L c we expect Ao > and hence for <p — to be 
locally stable. When L > L c , we have Ao < 0, and so 
ip = becomes unstable. In this case the ground state is 
given by a solution other than cp = 0. 

Suppose that we take the length scale L to be slightly 
greater than the critical value : 



L=(l + e)L c 



(3.2) 



with e > treated as small. We can scale the coordi- 
nates with appropriate factors of L to enable us to use 
dimensionless coordinates. The Laplacian —V 2 can then 
be expanded in powers of e using Eq. (3.2) to give 



- V 2 = -V 2 - eVi - e 2 V 2 



(3.3) 



where V 2 denotes that the Laplacian is evaluated with 
L = L Cl and Vi, V2, ■ ■ ■ are differential operators whose 
form is determined by the Laplacian and the specific cav- 
ity. (We will apply this to the rectangular box and the 
spherical cavity later to show how this works in practice.) 

An argument presented in Ref. Q, which is substanti- 
ated by application to specific exa mpl es, may be repeated 
to show that the solution to Eq. ( |2.5| ) can be written as 

ip(x) = e 1/2 [ipa(x)+e<px(x)+e 2 (p 2 (x)+---] , (3.4) 

where (po(x), ip \(x ), . . . are independent of e. Using 
Eqs. (|3.3|) and (3.4) in Eq. (2.1), with the potential given 



by Eq. (|j), and then equating equal powers of e to zero, 
we obtain a set of coupled differential equations, the first 
two of which are 











(3.5) 



Vi<A) - • (3.6) 



(Higher order equations may be obtained in a straight- 
forward manner if needed.) The set of differential equa- 
tions obtained may be solved iteratively be ginnin g with 
Eq. (f£j|. By comparing Eq. (f^) with Eq (]2 10|) where 
£q is defined by the critical length in Eq. (3.1), it is ob- 
served that ipo is the eigenfunction of the Laplacian whose 
eigenvalue is th e sm allest (with L set equal to L c ). 

Because Eq. (3.5) is a homogeneous equation, the over- 
all s cale of the solution is not determined. However 
Eq. (|3.6|), as well as the higher order equations not writ- 
ten down explicitly, are inhomo gen eous . The scale of ipa 
can be fixed by requiring Eq. (3.4) to minimize the en- 
ergy to lowest or der in e. F rom Eq. (2.6), using Eq. (3.4) 
along with Eqs. ( |3.5| ) and ( |3.6| ), it follows that 



E = I da T i, — a 
24 



where we have kept only the lowest order terms in e. If 



we let tpo be any solution to Eq. (3.S) which satisfies the 



correct boundary conditions, and write 

<po = Aifo 
for some constant A, we obtain 



E = C - C X A 2 + C 2 A 4 



where 



Co = J^ dax M a 
is independent of A and e, and 
1 



(3.8) 



(3.9) 



(3.10) 



Ci = |e 2 y der^oVi^o , (3-H) 
C2 = y/ £ da ^o , (3-12) 



are independent of A. For C\ > 0, we have 



A = ± 



Ci_ 

2C 2 



1/2 



(3.13) 



as the value of A which makes the energy a local mini- 
mum. This sets the scale of the solution for ifo : 



1/2 



with <po any solution to Eq. (|3.5[). 



(3.14) 
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To pr oceed to the next or der in e we must solve 
Eq. (3.6) with if given by Eq. ( |3.14 ). The general solu- 
tion may be expressed as 



(3.15) 



wher e ip ih is a solution of the homogeneous equ atio n 
Eq. (3.5) and tpi p is any particular solution to Eq. (3.6). 
The overall scale of <pi p is fixed since it satisfies an in- 
homogeneous equation; however, the scale of ipih is not 
determined. We can choose (pih to be proportional to 
ipo and fix the constant of proportionality by requiring 
that the energy be minimized to the next order in e. We 
therefore take 



= e 1 / 2 



Ci 



1/2 



— ) <p + e(A<p + ipip) 



+e 2 ip 2 



(3.16) 



in Eq. ( |2.6| ) and minimize the resulting expression for 
E(A). Because we are only interested in the value of A 
which makes E(A) a minimum, it is simpler to evaluate 

^(A) = e 3 / 2 £da^ (-VV+^ 3 -^V) (3.17) 

to lowest order in e. Because we have only solved the 
equation of motion up to and including terms of order 
£ 3/2 m obtaining Eqs. (3.5) and (3.6), the next term in 



the integrand of Eq. ( 3.17 ) will be of order e 5 / 2 . A short 
calculation shows that 



±E(A)=e\D a + D 1 A) 



where 
Do 

D 1 



2C 2 



1/2 



A fCi_ 
2 ^2C* 2 



(3.18) 
(3.19) 



'Po^Wo 



2C 2 



~4 
<fo 



We conclude that if D 1 > 0, then 
A 



Di 



(3.20) 



(3.21) 



gives the valu e whi ch minimizes E(A) to lowest order in 
e. From Eq. (pl| , in summary, we therefore have the 
approximate solution given to order e 3 / 2 by 



ip{x) 



M 2 



Cr\ 1/2 . D . 
2C 2 ^ + e(--^o- 



(3.22) 



where <p>o and <p\ p are any solutions to Eqs. (3.5) and 
( |3.6[ ) resp ectively. C\ and C 2 are defined by Eqs. ( 3.11) 
and (3.12). D and Di are defined by Eqs. (3.19) and 

( poD - 



It should be clear from the analysis that we have pre- 
sented how the method extends to any order in e. First 



solve the e qua tion of motion Eq. (2.3) to order e" +1//2 



using Eq. ( |3.4D extended to order ^+ 1 / 2 and Eq. (|3.3| ) 
to order e n . (n = 0, 1, 2, . . . here.) Evaluate -^E(A) to 



rder 



,2n+2 



using the solutions found, and then solve for 



A as we have illustrated for the cases n = 0,1. 



IV. THE RECTANGULAR CAVITY 
A. Exact result 

The simplest case of a cavity where an exact solution 
can be found occurs for a field confined inside of a rect- 
angular box. The first case we will look at is for a field 
which satisfies a Dirichlet boundary condition on oppo- 
site sides of one pair of box walls. We will choose this 
to be the y direction and take — L/2 < y < L/2 with 
ip = when y = ±L/2. In the x and z directions we will 
choose either periodic boundary conditions, or else Neu- 
mann boundary conditions with J^ip and J^ip vanishing 
at x — ±L x /2 and z = ±L z /2 respectively. With either 
of these two choices, the ground state will not depend on 
the x or z coordinates, and the problem reduces to one 
that is one-dimensional. One physical application of this 
is to the case of two parallel plates, as in the Casimir 
effect, where the plate separation is much less than their 
linear extent, (ie Keep L finite, and let L x , L z — > oo.) In 
this case, with L x , L z — > oo, the choice of boundary con- 
ditions in the x and z directions would not be expected to 
be important. Later in this section we will discuss what 
happens in the case where Dirichlet boundary conditions 
are imposed in all three spatial directions. 



If we assume ip — <p(y) only, then Eq. (2.5) becomes 



dy 2 6 



a 2 ) = . 



This equation admits a first integral, 



d,p 
dy 



2A W 



a'f - C , 



(4.1) 



(4.2) 



with C a constant. Note that ip — ±a solves Eq. flt.l|), 
but does not satisfy the requirement that the field va nish 
at (p = ±L/2, so is not allowed, p = satisfies Eq. ( |4.l| ) 
as well as the boundary conditions, so is a valid solution. 
To s ee if <p — is locally stable we look for the solution to 
Eq. ( 2.10 ) of lowest eigenvalue. Because ipo(y = ±L/2) = 
0, w e hav e tp oc sin (^{y + L/2)) and l\ = tt 2 /L 2 . From 
Eq. ( 2.11 ) we see that <p = is stable if L < L c where 



L r = 



1/2 



(4.3) 



If L > L c , then ip = is unstable, and because the 
boundary conditions prohibit constant values of ip, the 



5 



ground state must be spatially dependent. We will now 
find this solution. 

For tp(y) to be continuous on the interval [-L/2, L/2] 
with tp(±L/2) — 0, it must have a stationary value some- 
where. This allows us to conclude that C defined in 



Eq. (4.2) must be non-negative, and also that C < ^a 4 - 



We will define 



24 



(4.4) 



with w real and satisfying < w < 1. The cases w = 
and w = 1 require special treatment. Fo r w = 1 it is easy 
to show that the only solution to Eq. (4.2) which satis- 
fies the boundary conditions that p = at y = ±L/2 is 
if = for all y. We already know that this solution is un- 
stable for L > L c so we will concentrate on w < 1. The 
case w = is simple to solve, and leads to the usual kink 
solution. (See Ref. j|, [| for example.) This does not 
satisfy our boundary conditions. Therefore we will re- 
strict < w < 1. It is worth remarking on the difference 
between the confined case and the field in an unbounded 
region. In the latter case, the requirement that the field 
configuration leads to a finite energy requires the field 
to asymptotically approach a zero of the potential || |) . 
For a field confined by a cavity, any well-behaved func- 
tion will have a finite energy as a trivial consequence of 
the finite volume of the region; thus, the finite energy 
requirement does not give anything useful in our ca se. 

If we restrict < w < 1, then the solution to Eq. ( |4.2| ) 
can be expressed in terms of Jacobi elliptic functions M. 
We find 



<p(y) = ±a(l - W y/hn U (y + , (j-^j 



1/2N 



where 



= 



"I 1/2 



12 



(1 + w)a 2 



as the solution to Eq. ( |4.2| ) which vanishes at y = 
If we impose the other boundary condition at y 
we find the condition (ic^l) 



sn (3L, 



1 



1/2N 



= 



(4.7) 



Making use of the property that sn(u, k) has zeros 
for u = 2nK(k) for n = 0, ±1, ±2, . . . where K(k) is the 
complete elliptic integral of the first kind, 



K(k) 



we conclude that 



tt/2 



d0(l - k 2 sin 2 (j))- 1 ' 2 , (4. 



0L = 2nK 



1 — w 
1 + w 



1/2N 



(4.9) 



where n — 1,2,.... We can restrict n to non- negative 
values since sn(u, k) is an odd function of u, and a sign 
change of u just changes the overall sign of <p which is 
arbitrary. We rule out n = b ecause this corresponds 
to ip identically zero in Eq. (4.5), and we already know 
that this is unstable for L > L c . The interpretation of n 
is that (n — 1) gives the number of zeros of tp in the open 
interval {-L/2, L/2). 

The next question concerns the sol utio n for the arbi- 
trary integration constant w in Eq. (4.9). We need to 
solve 

7^d + -r 1/2 ^((^) 1/2 ) , (4.10) 



where 



La ( A 



ln ~ 2n I 12 



1/2 



(4.11) 



for w, with n = 1,2,.... The properties of the complete 
elliptic integral of the first kin d ma y be used to show 
that the right hand side of Eq. ( 4.10 ) is a monotonically 
decreasing function of w on the interval < w < 1, 
approaching infinity as w — > and the value _7r\/2/4 as 
w — > 1. This means that a solution to Eq. ( |4 . 1 0[ ) only 
exists if 



7n > 



ttV2 



(4.12) 



Making use of Eq. ( 4.11 ) and the definition of the criti- 
cal length in Eq. (4.3) shows that we will always have a 
solution for w in Eq. (4.10) if 



L > nL r 



(4.13) 



Apart from the different boundary conditions, the situa- 
tion is very similar to the twisted field case on the circle 
examined by Avis and Isham In fact, we can make 
use of their clever proof that for n > 2 the solutions for 
tp(y) are all unstable. This is not unexpected, since the 
energy of solutions with n > 2 are all greater than that 
for the n = 1 solution. We will therefore concentrate on 
the case of n = 1. 

If we u se the definition of L c given in Eq. (4.3) in 
Eq. ( 4.10| ), when n = 1 we find 



nL 



2^2L 



Vl + w = K 




(4.14) 



showing that the solution for w depends only on the di- 
mensionless ratio of L/L c . By expanding both sides of 
Eq. (4.14) in powers of w, it is easy to show that for large 
values of L/L c a good approximation for the solution is 
given by 



w ~ 8 exp 



irL 
V2L C 



(4.15) 
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Even for relatively small values of L/L c this turns out to 
be quite accurate. (For example, when L/L c = 3, we find 
an agreement to six decimal places b etwee n this approx- 
imation and the resu lt of s olving Eq. ( 4.14 ) numerically.) 
The solution to Eq. (4.14) is shown in Fig. 1 for a range 
of L/L c . It is clear from this figure that as the value of 




L/L 



FIG. 1: The solution for w obtained by solving Eq. (4.14) as 
a function of L/L c . 

L/L c is increased, w tends towards 0, as predicted from 
the approximation Eq. (4.15). Since w = corresponds 
to the constant solution ip = a, it would be expected that 
as we increase the value of the ratio L/L c , the solution 
for ip will try to be as close as it can to the constant value 
of a. This is in fact what happens as we show in Fig. 2. 
As L/L c gets larger the solution tends towards the step 







/ (p/a 


^^""--^ L/L c =$\ 


/ / 08 


L/L=2\ \ 


/ / 06 


zyz^i.i \ \ 


/ / 0,4 




/ / / 0.2 


L/L =1.01 \ \\ 




L/L = 



FIG. 2: The solution for <p(£)/a is shown for various values 
of L/L c as a function of £ = 2y/L. 



function, which is as close as the boundary conditions 
allow it to get to the standard solution ip = a. 



B. Approximate result 

As a test of the approximation method described in 
Sec. Ill, we will use it to compare with the exact solu- 



tion we have just fou nd. Take L to be close to the crit ical 
length L c as in Eq. (3.2) with L c defined by Eq. (4.3) in 
the present case. As explained in Sec. H, it is advan- 
tageous to adopt dimcnsionless coordinates, so we will 
define 



(4.16) 



with -1<£<1. Since = — -j, by using Eq. (4.1C) 



dy 

and the expansion of L about L c as in Eq . ( |3 . 2| ) , it is easy 
to read off Vi and V2 defined in Eq. (3.3) to be 



V 2 = — 



12d?_ 



We also have 



V, = 



±^1 

li de 



(4.17) 

(4.18) 
(4.19) 

(4.20) 



The next step is to solve Eq. (3.5) for any solut ion (pp. 
Using Eq. ( 4.2C ) and the definition of L c in Eq. (P|), it 
is easy to see that 



Po(0 = cos (|^) 



(4.21) 



is a solution which satisfies the proper boundary condi- 
tions. A straig htfor ward calculation of C\ and C2 defined 
by Eqs. ( |3 . 1 1 ) and Q3.12 ) leads to the conclusion that 



2C 2 



r 



(4.22) 



We therefore have as our leading order approximation to 
the ground state 



p(0 



3 6 



1/2 



a cos 



(4.23) 



Before comparing this with the expansion of the exact 
resul t, we will first evaluate the next order c orre ction to 
Eq. (|4.23D , This entails initially solving Eq. (Q for any 
solution pip. With 



MO = ( 3 



1/2 



a cos 



(4.24) 



and Vi given by Eq. ( 4.17 ), it can be shown that a par- 
ticular solution to Eq. (3.E) is given by 



(4.25) 
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The constants Dq and D\ defined in Eqs. (3.19) and 
( 3.20 ) may now be evaluated with the result 



Do 



17\/6 
-36- 



(4.26) 



Combining all of these results, and using Eq. ( |3.22 ), we 
obtain the approximate ground state solution to be 




3tt 

a cos I — f 
18 V 2 ^ 



, (4-27) 



up to, and including, terms of order e 3 / 2 . (We drop the 
± here.) 

The e xact ground state solution was found to be given 
by Eqs. (|4.5|j4.6| , H) with n = 1. Using the dimensionless 
coordinate £ defined in Eq. ( 4.16| ) the exact solution reads 



p(0 = o(l - W ) 1 / 2 sn((l + 0K(k w ), kj) , (4.28) 
where 



1 — w 
1 + w 



1/2 



(4.29) 



We will now expand this result consistently to order e 3 / 2 
with L = (1 + e)L c . 

For L close to L c , w will be close to 1; thus we will let 

w = l-r), (4.30) 

with rj << 1. T he fir st tas k is t o determine rj in terms 
of e. From Eqs. ( 4.10 ) and ( 4.11 ), with n = 1, it can be 
shown that 



1 



ttV2 



(i + w )-^ 2 K(k w ) 



(4.31) 



(The result has been simplified using Eq. (4.3).) Using 
Eq. (4.30) in Eq. ( 4.29Q , and expanding the complete 
elliptic integral of the first kind in powers of rj results 
in 



1 



1 



57 
256 



This expansion may be inverted to find 



n 



38 



(4.32) 



(4.33) 
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which when used in Eq. ( 4.30] ) , gives us an expansion of 
the parameter w in terms of e. 

The final task is to expand the Jacobi elliptic function 
in powers of k w when k w is small. A useful result is 
contained in Ref . || , and reads 

sn(u, k w ) ~ sin it + -e(sin it cos u — u) cos it , (4.34) 



in our case. The final step is to note that from Eq. (4.2? 
we have 



(l + 0K(k w ) 



(4.35) 



so that u in Eq. (4.34) also depends on e. If we expand 
to order e 3 / using the expansions just described, 
after a short cal culat ion we obtain a result in exact agree- 
ment with Eq. ( [4.271) . 

The main conclusion of this section is that the ap- 
proximation method is in complete agreement with the 
expansion of the exact result, at least for the first two or- 
ders in the expansion used. This is sufficient to generate 
some faith in the general procedure outlined in Sec. Ill in 



cases where it is not possible to find an exact solution by 
analytical means. An example will be given in the next 
subsection. It can also be used to provide a useful check 
on the results of numerical calculations. We will study 
such a case in Sec. ^. 



Dirichlet boundary conditions in three 
dimensions 



In this section we will examine the case of a cubical 
box of side length L with the field vanishing on all of 
the box walls. Although this is perhaps a more realis- 
tic situation for a confined field than that considered in 



Sec. IV A, unfortunately we have not been able to find 
the exact solution when the zero field is unstable. This 
is quite unlike the situation where periodic or antiperi- 
odic boundary conditions are imposed on the walls. In 
this case the exact solution in three dimensions can be 
simply related to that found in one dimension ||. The 
Dirichlet boundary conditions thwart the application of 
a similar procedure here. However we can still use the 



approximation method described in Sec. Ill 



The stability of the solution ip = is determined by 
the lowest eigenvalue £q in Eq. (2.10). It is easy to show 
that 



, , . /irx\ /ttv 

ipo(x,y, z) = cos J cos y— | cos 



(?) («•) 



is the eigenfunction of the Laplacian with the lowest 
eigenvalue given by 



£ 0- L2 



(4.37) 



The critical l ength L c is the value of L for which Ao de- 
fined in Eq. (2.11) vanishes. This gives 



L = 



18tt 2 
Aa 2 



(4.38) 



For L < L c , ip = is stable, while for L > L c it is unsta- 
ble. Thus when L > L c the ground state will involve a 
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non-constant value of ip in order to satisfy the Dirichlet 
boundary conditions. 

It proves con venie nt to define dimensionless coordi- 
nates as in Eq. (4.16) : 



y : 



2 6 ' Z 



§6 



(4.39) 



so that the bound ary of the cub e is at £j = ±1 for i 
1, 2,3. From Eqs. $3) and fO) we find 



Vi 



V, = 



± V — 

H V — 



(4.40) 
(4.41) 
(4.42) 



A solution to Eq. (3.5) with the correct boundary condi- 
tions is 



^o(£i,6,&) = n c ° s (f & 



i=l 



(4.43) 



This ma y be used to calculate C\ and Ci in Eq. (3.11) 
and Eq. ( 3.12| ) resulting in 



Ci_ 

2C 2 



128 

-w 1 



(4.44) 



The leading order appro ximation to the ground state 
therefore, from Eq. (3.14), involves 



<£>o(£i>62>&) 



8^6 



a II C0S (f&) 



(4.45) 



To proceed to the next ord er we mus t solv e the partial 
differential equation Eq. ( |3.6| ) with Eq. (4.45) substituted 
for ifQ. A solution with the correct boundary conditions 
can be shown to be 



<Pip 



2^6 
~8T 



a {r os (T a ) cos (T 6 ) cos (T^) 



+9 



2 6 J cos It 6 J cos lT ej 

+ (1 <-► 2) + (1 <-> 3) 
+ (1 ^ 3) + (2 ^> 3)]} 



3tt 



(4.46) 



where, to save space, we have used (i <-> j) to mean the 
first term in the square brackets with the indices i and j 
on £j and £j switched. 

Finally we evaluate Do and Di defined by Eq. ( 3.19| ) 
and Eq. (3.20). A straightforward calculation leads to 



D Q _ 1375V6 
ITi ~ 4374 



(4.47) 



The solut ion may now be wr itten down immediately from 
Eq. (3.22) using Eqs. flOg) and j07| ) with (fo given n 
Eq. ( 4.45 ) and tp± p in Eq. ( [4.46 ). (To save space we will 



not write this out explicitly.) 



V. THE SPHERICAL CAVITY 

We now examine the case where the scalar field is con- 
fined by a spherical shell of radius R, with the field van- 
ishing at r = R. The ground state should be spherically 
symmetric, so that ip will only depend on the radial coor- 
dinate r if we use the usual spherical polar coordinates. 
We need to solve 



rv2 A 



where (p = ip(r) and 



6 



1 d 

r 2 dr 



2 2\ 

a 



0, 



dr 



(5.1) 



(5.2) 



with cp(r = R) = 0. 

tp = is obviously a valid solution, and from the discus- 
sion in Sec. II, the stability is determined once we know 
the lowest eigenvalue of the Laplacian as in Eq. (2.10). 
In the present situation it is easy to show that 



-00 



sin (nr/R) 



is the eigenfunction of lowest eigenvalue, with 



f- 



IT \ - 



(5.3) 



(5.4) 



From Eq. (2.11), we can conclude that ip — is stable for 
R < R c where the critical shell radius R c is defined by 



c \ 1/2 
6 \ 7T 



(5.5) 



For R > R c , <p = is not the ground state. 

To find the stable ground state when R > R c , we must 
solve the non-linear differential equation Eq. (5.1). We 
were unable to find an exact analytical solution here, un- 
like the case of the one-dimensional box where the dif- 
ferential equation was able to be solved by quadrature. 
Instead we solved the equation by numerical integration. 
The results are plotted in Fig. 3 for a range of values 
of R/R c . The results show that as the value of R/R c 
is increased, the value of the field tries to get closer and 
closer to the constant value of ip = a over a greater range 
of values of the radius. Because the boundary conditions 
require <p(r = R) = 0, this results in a very sharp drop- 
off in the value of the field as r — > R. The profile of the 
field starts to approach a step-function as R/R c — * oo. 

Although we were not able to find an analytic solution 
to Eq. (5.1), we can still apply the approximation method 
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e 


<p(f})/a 


order e 1//2 


order e 3 ' 2 


0.001 


0.068291 


0.068370 


0.068288 


0.01 


0.213584 


0.216206 


0.213590 


0.1 


0.608397 


0.683703 


0.601004 


0.2 


0.774283 


0.966902 


0.732994 


0.3 


0.862136 


1.184208 


0.754491 



TABLE I: A table showing the value of the field at the o rigi n 
in units of a. e shows how close R is to R c as in Eq. (5.7). 



The second column gives the true value of <p(0)/a found by 
numerical integration of Eq. (5.1). The final two columns 
show respec tivel y the results found from the approximate so- 
lution Eq. ( 3.22 ) using the leading order term, and the next 
order correction. 



FIG. 3: The solution for <p(r)/a is shown for various values 
of R/R c as a function of r/R. 



to be used in Eq. (3.6). 

A particular solution to Eq. ( |4.6| ) can be shown to be H 



described in Sec. III. A dimensionless radial coordinate 
p is denned by 



and we set 



r = Rp, 



R=(l + e)R c 



(5.6) 



(5.7) 



as in Eq. (jTg). With V 2 given by Eq. (|j) we find the 
expansion Eq. (|3.3|) where 



V 2 . 
Vi 

v 2 



— V 2 



v , 



Rl o 

3 rv2 

— V 

m p ' 



(5.8) 
(5.9) 
(5.10) 



and V 2 given by taking r = p in Eq. (5.2). 

The lowest order contribution to the approximate so- 
lution follows from solving Eq. (ETfl). This results in 



<Po{p) 



sin(7rp) 



(5.11) 



The overa ll scal e is set by calculating C\ and C2 in 
Eqs. (|3.1l|) and (|3.12|) with the result that 



2C 2 



1/2 



a 



with |] 



This gives 



J=-[Si(27r)-Si(47r)] 



<Po(p) 



a sin(7rp) 

VI p 



(5.12) 



(5.13) 



(5.14) 



Vip(p) 



r [Ci(4 7 rp)-Ci(2 7 r j0 )] 



sin(7rp) 

2/3/2 l^h"*»av - ^\«»i*n — - — 
r [Si(47rp)-2Si(27rp)]- COsM 



2J 3 / 2 ' 

■J1J2 cos w) 



P 



(5.15) 



The compli cated nature of this renders the calcul ation of 
D in Eq. fl3.19| ) difficult, although D x in Eq. fl3.20| ) is 
easily evaluated. It can be shown that 



D 

~ 0.417 349 770 a 

Di 



(5.16) 



This is sufficient to determine the approximate solution 
in Eq. ( |3.22[ ) correct to order e 3/2 . 

As a check on the numerical results shown in Fig. 3, we 
plotted the approximate solution we have just described. 
For small e the result was found to be indisti ngui shable 
from the result of numerical integration of Eq. (5.1); how- 
ever, as e is increased the agreement becomes less good 
as would be expected for a small e expansion. The dis- 
agreement is largest near p — where the field has its 
largest value. To get an idea of how close the approx- 
imate solution is to the true result, some of the values 
found for the field at the origin are included in Table 1. 
It is clear that by including the correction to the leading 
term a more accurate result is obtained for a wider range 
of e. Even for relatively large values of e (e ~ 1) where 
the small e expansion would not be expected to be par- 
ticularly accurate, we found the approximate solution to 
be close to the true one for p close to p — 1 , correspond- 
ing to the radius of the shell. Thus the approximation 
method described in Sec. ^TJ can provide a useful check 
on the results of numerical integration in cases where no 
exact solution is known. 
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VI. SUMMARY AND CONCLUSIONS 

We have presented an analysis of the ground state for a 
real scalar field confined in a cavity. Although we concen- 
trated on the double- well potential (2.2) it would be easy 
to extend the analysis to other potentials. (For example, 
in the case the bi-cubic potential, an exact solution in the 
cavity can be found in terms of the Weierstrass elliptic 
function as in Ref. In addition, the analysis can be 

extended to curved space, although it would be difficult 
to find exact solutions except in special cases. 

We presented a method for obtaining approximate an- 
alytical solutions in the case where the length scales of 
the cavity were close to the critical values at which ip = 
became unstable. For the rectangular cavity, we showed 
that this approximation method agreed with an expan- 
sion of the exact solution which we found. When ip = 
was not the ground state, it was shown that as the ratio 



of the size of the cavity to the critical size L/L c increases, 
the ground state tends towards the constant value which 
minimises the potential over as much of the cavity as pos- 
sible. The condition that the field vanishes on the box 
walls means that the field must always drop to zero, with 
the drop-off becoming increasingly sharp as L/ L c is in- 
creased. A similar behaviour was found for the case of a 
spherical cavity. By extrapolation, for a field confined to 
vanish on the walls of a general cavity, if the cavity is suf- 
ficiently large it would be expected that the ground state 
would correspond to a field which was constant almost 
everywhere inside the cavity, with a very sharp drop-off 
to zero as the cavity boundary is approached. 

Although we have concentrated on fields confined by 
cavities in the present paper, the general method de- 
scribed above can be useful in other situations where the 
boundary conditions prohibit constant values of the field. 
We plan to report on this elsewhere. 
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